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ABSTRACT 

We measure the projected rotational velocities of the late B-type hyperveloc- 
ity stars HVS7 and HVS8 from high resolution spectroscopy to be 60 ± 17 km s" 1 
and 260 ± 70 km s -1 . The 'slow' rotation of HVS7 is in principle consistent with 
having originated in a binary system, assuming a high inclination angle of the 
stellar rotation axis. However, the fast rotation of HVS8 is more typical of single 
B-type stars. HVS8 could have therefore been ejected by a mechanism other 
than that proposed by Hills. We also estimate the effective temperatures and 
surface gravities for HVS7 and HVS8 and obtain an additional measurement of 
their radial velocities. We find evidence in support of a blue horizontal branch 
nature for HVS7, and a main sequence nature for HVS8. 

Subject headings: stars: early-type — stars: fundamental parameters — stars: 
rotation — Galaxy: halo — Galaxy: stellar content — Galaxy: center 



1. Introduction 



The recent discovery of 10 hypervelocity st ars (H VSs; iBrown et al.ll2005l : lEdelmann et al. 
20051 Ifflrsch et allbood : brown et al]l2006al lbl. [ 20071 ) has raised many question s abo u t thei r 
nature and origin. The most widely accepted ejection mechanism, proposed by iHillsl (119881 ). 
involves the encounter of a close binary with a supermassive black hole (SMBH). Other 
possible mechanis ms ejecting stars from the Galactic center involve interm ediate- mass black 



holes (IMBHs; e.g. Yu fc Tremainel 



2003 



black hole (BMBH; e.g. lYu fc Tremainel 120031 : Merritt 



2006 



Lockmann fe Baumgardtl 120081). a bina r y mas sive 



Sesana et all 12008 . 120071 ). or 
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a clus ter of stellar mass black holes around the SMBH (jO'Leary fc Loebl 120081 ) . [Hansen 
( 120071 ) claimed that the rotational velocities of HVSs should be lower than those measured 
for single stars of the same spectral type if they originated in binaries, because of tidal 
effects. He predicted that the rotational velocitie s of the known B- t ype H VSs should be 
~ 70 — 90 km s" 1 , based on values com piled by lAbt fc Boonyaraki (12004 ) for B-stars in 
binaries. lLockmann fc Baumgardtl (120081 ) predicted high rotational velocities for HVSs that 
were ejected by a very close encounter with an IMBH in the Galactic center, however such 
encounters are very unlikely. 

These predictions cannot be tested with existing observations, as the low resolution of 
the discovery spectra of most HVSs is not sufficient to determine projected rotational veloci- 
ties (vsini). The onl y HVS with high resolut ion spectroscopy and a vsini mea surement is HE 



20081). in agreement with the prediction of 



0437- 5439, found by lEdelmann et al.l (12005). It ha s vsini = 55 ± 1 k m s 1 (iBonanos et al. 



Hansen ( 2007 ). However, Bonanos et al. ( 20081 ) 



and iPrzybilla et al.l (120081 ) measured half-solar metallicity for this early B-star, establishing 
its origin in the Large Magellanic Cloud (LMC). The possible ejection mechanisms for this 
star include an interaction with an IMBH or a SMBH, and a dynamical interaction of a 
single star in a dense cluster. This example demonstrates the importance of high resolution 
spectroscopy for understanding this newly discovered class of objects. 



Of the remaining HVSs, HVS2 (or US 708; iHirsch et al.ll2005l ). is classified as an evolved 
sdO star and reasonably well understood. However there is some ambiguity in the nature of 
the late B-type HVSs, since at their temperatures and gravities, the blue horizontal branch 
(BHB) crosses the main sequence. Hot BHB s tars generally have low rotational velocities 
and peculiar chemical abundances (jBehrll2003al ). thus high resolution spectroscopy of these 
faint HVSs can determine their nature by measuring their atmospheric parameters, chemical 
abundances and vsini. In addition, time series photom etry can reveal pulsa tions and confirm 
their main sequence nature, as was done for HVS1 by iFuentes et al.l (120061 ). 



Motivated by the lack of vsini and stellar parameter measurements for most of the 
known HVSs and the possibility of testing the nature of the SMBH in the center of our 
Galaxy, we performed high resolution spectroscopy of two HVSs. In this Letter we present 
our results. 



Observations 



We collected spectra of HVS7 and HVS8 (S PSS J113312.12+0108 24.9 and J094214.04+200322.1) 
with the blue chip of the MIKE spectrograph (jBernstein et al. 2003 ) installed at the 6.5-m 
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Magellan Clay Telescope at Las Campanas Observatory (Chile), on two half nights on UT 
2008 January 18-19. Each star was observed twice, with individual exposure times between 
900 and 1200 s, using a 1" x 5" slit and 3x3 binning. The total exposure times were 2100 s 
for HVS7 and 2400 s for HVS8. The resolution of the spectra is R = 32,000 at 4500 A. 



Th e spe ctra were extracted using the MIKE reduction pipeline (IKelson et al.l 12000 



Kelson! 120031 ) . The extracted spectra for each star were then averaged, normalized and 
merged. The wavelength coverage of the merged spectra is 3900-5050 A, with an average 
S/N-ratio per pixel of 15 for HVS7 and 14 for HVS8, based on the extracted continuum 
around 4500 A. These S/N-ratios and our spectra l resolution are sufficient to distinguish 



between high (> 130 — 140 km s ; lAbt et al 



2002) and low (70 - 90 km s _1 ; lHansenll2007h 



vsini values for B-stars. Next, we corrected the wavelength scale for Doppler shift, to allow 
comparison of the spectra with models (see §3). We measured t he heliocentric radial velocity 
of each star using the IRAF0 cross-correlation package RVSAO (IKurtz fc Minklll998l ) and the 
grid of models desc ribed in §3. Table [TJlists our results and the values previously reported by 
Brown et al.l (120071 ). §4 discusses the implications of our new radial velocity measurements. 



3. Spectral Analysis 



Our high resolution spectra allow direct determination of the effective temperature T e ff, 
surface gravity logg, and vsini of the stars by comparing synthetic model spectra to the 
observations. The S/N-ratio of the data is however too low to reliably measure abundances. 

We gene rated a g rid of synthetic spectra using the LTE ATLAS9 models and opacities 
developed by lKuruca (119931 ) . The grid covers T e ff between 8000-15000 K in steps of 1000 K, 
and logg between 3.0-5.0 in steps of 0.25 dex. The metallicity was set to solar, assuming 
that the HVSs are ejected from the Galactic center, where abundances are solar or supersolar 



(ICarr et al.l l2000l : iRarmrez et al.l I200CH : iNaiarro et al.l 120041 ; IWang et al.l 120061 : ICunha et al. 



20071 ). For the macro- and micro-turbulence velocitie s we adopted and 2 km s -1 , which 
are typical for late B-stars (IFitzpatrick fc Massalll999l ). The models were broadened by 0.15 
A to match MIKE's instrumental profile and resampled to a dispersion of 0.03 A/pix to 
match the dispersion of the stellar spectra. Finally, we convolved each model with rotational 
profiles between 10-350 km s^ 1 in 10 km s _1 velocity increments. 



Simultaneous fits to T e *f, logg and vsini were performed for each star by iteratively 



1 IRAF is distributed by the NO AO, which is operated by the Association of Universities for Research in 
Astronomy Inc., under cooperative agreement with the NSF. 
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comparing each model to the data. The agreement between each model and the observed 
sp ectra i s quant ified by the spectroscopic quality-of-fit parameter, z (normalized x 2 ), defined 
by iBehrl (l2003al ) and given by the equation 



A r , 



points 



rms 2 



2 

rms min 



where N po i nts is the number of points in the spectrum, and rms and rms m i n are the root 
mean squared deviation between each model and the stellar spectrum, and the smallest 
value of the rms found, z = gives the best model fit, and z = 1 defines the statistical la 
confidence interval of the result. The following subsections describe the derivation of T e ff, 
logg and vsini for each target. 



3.1. HVS7 

The spectrum of HVS7 (V=17.80 mag, S/N=15) includes four Balmer lines {Hp - H e ) 
from which T e ff, logg and vsini can be estimated. We also detect Mg II, Si II, Ca II and 
Fe II lines that can in principle be used to further constrain T e ff and vsini, however the two 
main T e ff indicators (Mg II and Si II) have anomalous line strengths and cannot be used to 
constrain the T e ff. We are therefore left with only the Balmer lines that are simultaneously 
sensitive to T e ff, logg and vsini, bu t can still provide non -degenerate values of these param- 



eters for late B-type stars (see §4 of lMcSwain et al.ll2007l ). We performed two tests to verify 



that Balmer lines alone are sufficient to simultaneously derive the three parameters: a) we 
applied our analysis to synthetic spectra with added random noise matching the S /N-ratios 
of the observations, and in all cases recovered the input values within errors, b) we applied 
our analysis to a high S/N spectrum of a late B-type star (HR7447, B5 III), kindly provided 
by L. Lyubimkov. Our analysis yielded T e ff = 14,000 ± 1000 K, logg = 3.75 ± 0.25 dex, 



and v s ini = 70 ± 20 km s , in agreement with the parameters derived by lLyubimkov et al. 



(120021 . 12004! . T ef f = 13,400 K, logg = 3.64 dex, vsini = 76 km s" 1 ). 



We then proceeded to fit the spectrum of HVS7 for T e ff, logg and vsini. We performed 
several tests to determine the stability of the best fit solutions. We ran fits to the entire 
spectrum, 100 A windows centered on each Balmer line (to ensure that the wings and 
some continuum are included), portions of the spectrum outside the Balmer lines, and 10- 
20 A windows centered on metal lines. In the last two cases we had to fix T e jf and logg to 
the values from the fits to the entire spectrum and the Balmer lines and only fit for vsini. All 
the tests give fully consistent results, with the following best fit parameters: T e ff = 12,000 
± 1000 K, logg = 3.50 ± 0.25 dex, and vsini = 60 ± 17 km s _1 . The z minimization results 
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for the full spectrum are shown in Figured! We have adopted conservative errors for T e ff 
and logg equal to the grid step size, versus their smaller la statistical errors. The statistical 
la errors for vsini (horizontal dotted line in the z vs. vsini plot in Figured]), are —20 and 
+50 km s _1 , however visual comparison of the models to the observed spectrum show they 
are too large. Instead we adopted the errors resulting from the fits to individual metal lines. 
The left panel in Figure [2] compares the best fit model to the Balmer lines of HVS7. 

Figure [3] shows metal lines detected in HVS7 with vsini = 40, 60, and 80 km s _1 models 
overplotted. The purpose of this plot is two-fold; the left-side panels show how the Ca II 
K and Fe II A4233 lines (vsini is derived from these two and the Fe II AA4549 and 4583 
lines) give vsini = 60 km s" 1 as the best model fit, and their depths agree with the solar 
abundance adopted in the models. The vsini from these lines also agrees with the fits to 
the full spectrum and the Balmer lines. The right-side panels show the behavior of the Si II 
doublet, the Mg II and He I lines. None of the models in the grid can reproduce the depths 
of those lines. While He I and Mg II seem depleted in the atmosphere of HVS7, the Si II 
4128/4130 lines seem strongly enhanced. The models cannot reproduce either the depth nor 
the line strength ratio of the He I and Mg II lines. The enhancement of the Si II 4128/4130 
doublet is even more significant when taking into account that the Kurucz LTE ATLAS9 
models overpredict t he strengths of thes e lines. This problem persists even after including 
non-LTE corrections (jSmartt et al.ll200lf) . Abundance peculiarities have been noticed before 



in BHB stars by several authors (IGlaspey et al.lll989l : iMoeh 
however, most of these are very slow rotators (< 8 km s" 



-l. 



er et al 



19991 : iBehr et al.lll999l ): 



Behrll2003al ). 



3.2. HVS8 

The spectrum of HVS8 (V=18.09 mag, S/N = 14) has a S/N similar to the HVS7 
spectrum, however, inspection of the spectrum of HVS8 for metal lines gives null results. 
This can be explained by very low metal abundances, strong depletion, or highly broadened 
metal lines. The vsini obtained below points towards the latter case. 

To derive T e ff, logg and vsini we used only the spectrum above 4000 A because of 
problems with the continuum normalization at shorter wavelengths. As with HVS7, we 
simultaneously fit for T e ff, logg and vsini by iteratively comparing the spectrum of HVS8 
to our model grid. We ran fits to the entire spectrum and 160 A windows centered on the 
Balmer lines, which gave consistent parameters: T e ff = 11,000 ± 1000 K, logg = 3.75 ± 
0.25 dex, and vsini = 260 ± 70 km s _1 . The lack of metal lines in the spectrum of HVS8 is 
consistent with a high vsini that results in strong line broadening. The error in vsini in this 
case comes directly from the z — 1 statistical la result, as visual comparison of the spectrum 
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to the models does not allow us to place a finer constraint. The right-side panel in Figure [2] 
compares the best fit model to the Hp, H 7 and H$ lines of HVS8. The flat-bottomed cores of 
the Balmer lines, which are the most sensitive regions to vsini, clearly illustrate that HVS8 
rotates faster than HVS7. 



Radial Velocities 



The new radial velocity observations in Table [T] provide a third epoch for each star 
and allow to check for variations. Our radia l velo city measurement fo r HVS7 is identical to 
the values reported by iBrown et al.l (j2006bl ) and iBrown et al.l (120071 ). within errors. Such 
measurements provide clues to the nature of HVSs. As pointed out by Brown et al., de- 
termining the nature of late B-type HVSs is not straightforward because late-type main 
sequence B-stars and hot Blue Horizontal Branch (BHB) stars have identical atmospheric 
parameters. BHB stars are less luminous and therefore closer in distance to us. Establishing 
the evolutionary stage of HVS7 is critical because its radial yelocity is marginall y consistent 
with it being a BHB runaway star bound to our Galaxy (IBrown et al.l l2006bl ). The lack 
of significant radial velocity variations for HVS7 suggests it is not a binary, nor a pulsator. 
Slowly pulsating main s equence B-type stars typically show ra dial velocity variations of ~ 
20 km s _1 in amplitude (lAerts et al.lll999l ; iMathias et al.ll200lf). while BHB stars appear to 
be st able, as they fall outside the RR Lyrae instability strip (jContreras et al.ll2005l ; ICatelan 
20051 ). The long term radial velocity stability of HVS7, combined with the metal abundance 
anomalies (see §3.ip . hint towards HVS7 being a BHB star. Its vsini (60±17 km s _1 ) is 
higher than typically found for BHB s tars (< 8 km s _1 ), although rotators with vsini up to 
~ 40 km s _1 have been observed (e.g. Behr 2003al lbl). The true nature of HVS7 as a bound 
BHB star will have to be disentangled by astrometry. 

For HVS8 we detect a radial velocity variation of 23 km s" 1 , consistent with a pulsating 
main sequence B-type star. We cannot discard the possibility of HVS8 being a binary, 
although the system would have a very low mass-ratio, since there is no evidence of lines 
from a com panion in the spectr um. The star is most likely a main sequence slow pulsator, 
like HVS1 (IFuentes et al.l I2006T ) . The low S/N of our spectrum does not allow to test for 
metal abundance anomalies in HVS8, however, the high rotational velocity of this star will 
make its abundance analysis difficult, even with higher S/N spectra. 
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Discussion 



We have derived vsini, T e ff and logg for HVS7 and HVS8, two of the ten currenty 
known HVSs. Th eir T r jf and l o gg are consis tent with the stars being late B-type, as ini- 
tially classified by iBrown et al.l (j2006bl . 120071 ) using photometric color indexes. HVS7 has 
a projected rotational velocity vsini = 60 ± 17 km s -1 , while for HVS8 vsini = 260 ± 70 
km s -1 . These measurements provide the first direct observational test to the prediction by 
Hansen! (120071 ). who suggests that HVSs ejected via Hills' mechanism should rotate system- 
atically slower ( 70 — 90 km s _1 ) than single stars of the same spectral type (134 ± 7 km s _1 ; 
Abt et al.l 120021 ). If the HVSs have fast rotational velocities typical of single B-type stars in 
the field, other ejection mechanisms, such as three-body encounters of single stars with ~ 
10 3 -10 4 M IMBHs, with a binary MBH, or with - 10 M stellar -mass black holes orbiting 
the Galactic SMBH, have to be invoked. 

The vsini values of HVS7 and HVS8 are lower limits to their true rotational velocities, 
imposed by the inclination angle of the rotation axis of the stars. If the inclination of the 
rotation axis of HVS7 is low, its rotational velocity could in principle be much higher. In 
that case both targets are inconsistent with Hansen's prediction for Hills' scenario. However, 
a sample of only two vsini measurements is not enough to conclusively discern between the 
different scen arios proposed and more vsini measurements are necessary. Statistical tests 
performed by iPeretd (120071 ) conclude that a sample of 25 or more HVSs will be needed to 
distinguish between scenarios at a > 95% confidence level. 

We also detect abnormal enhancement and depletion effects in the strength of some of 
the metal lines of HVS7. These anomalies, together with the apparent lack of pulsations 
are consistent with HVS7 being a BHB star. Its vsini is also marginally consistent with it 
being a fast-rotating BHB star, if the inclination angle of the star's rotation axis is close to 
90 degrees. However, confirmation that HVS7 is a BHB star will not be possible until its 
proper motion is accurately measured. Finally, we find evidence of radial velocity variations 
in HVS8 consistent with a pulsating main-sequence B-type star nature. Additional radial 
velocity measurements and time-series precision photometry will confirm this detection. 
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Fig. 1. — z minimization result for full spectrum of HVS7. Left: Contour plot of logg vs. 
T e ff for a fixed vsini of 60 km s -1 ; contours correspond to z — 1, 2, 3, 4; crosses to the models 
in the Kurucz grid. Right: z vs. vsini for the best fit values T e ff = 12,000 K and logg = 
3.50 dex. z — 1 shows the statistical la error of the fit; the result is vsini = 60 ± 17 km s -1 . 
A similar analysis was done for HVS8. 
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Table 1. Heliocentric Radial Velocities 



Target 


RV (km s- 1 ) 


RV (km s- 1 ) 




(this work) 


(previous work* ) 


HVS7 


529 ± 2 


531 ± 12 


HVS8 


489 ± 2 


512 ± 10 



IBrown et al.l (120071 ) report two radial ve- 
locity observations for each star and find iden- 
tical values within their errorbars. 
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Fig. 2. — Comparison of the best fit models to the Hp, H 1 and H$ Balmer lines of HVS7 
(left) and HVS8 (right). Models are shown as red (dotted) lines. 
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Fig. 3. — Metal lines detected in HVS7; the overplotted models correspond to T e jf = 12,000 
K, logg = 3.50 dex and vsini = 40 km s _1 (blue), 60 km s" 1 (red), and 80 km s _1 (green). 
Ca II A3933 and Fe II A4233 are well reproduced by our best fit model (red), while Si II 
AA4128-30, He I A4471 and Mg II A4481 show abundance peculiarities. The sharp absorption 
line at 3926.5 A corresponds to interstellar Ca II K after correcting the radial velocity of 
the HVS7 spectrum to the heliocentric reference system. 



